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(57) The present invention provides a signal 
processing apparatus and method for enhancing the 
communication quality and increasing the communica- 
tion capacity by reducing the effects of interference and 
noises with the nice beam pattern. And, the inventive 
signal processing apparatus and method introduce a 
simplified computational technique for generating the 
nice beam pattern having its maximum gain along the 
direction of the wanted signal and maintaining the gain 
toward the direction of the interfering signals in as low a 
level as possible. 
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Description 

FIELD OF THE INVENTION 

This invention relates to a signal processing technique for wireless communication systems, and more particularly 
to a signal processing apparatus and method for reducing the effect of interference and noise by controlling beam pat- 
terns in real-time, in a communication system utilizing an antenna array. 

BACKGROUND OF THE INVENTION 

In general, an original signal transmitted by a certain transmitter (hereinafter, simply called a "wanted signal") is 
always received at a receiving set together with other plural interfering signals. 

Since the level of distortion in a telecommunication system is determined by the ratio between the power of the 
wanted signal and total power of all the interfering signals, even if the level of the wanted signal is much higher than 
each of the interfering signals, the distortion of the communication system can pose a serious problem when the total 
power of all the interfering signals proportionally increased according to the number of the interfering sjgnals is rather 
high. In conventional telecommunication systems, interfering signals make it very difficult to extract the information from 
the wanted signal. 

Although an antenna array system has been considered as a countermeasure to improve the problems caused by 
the interfering signals, no practical method of synthesizing the antenna array system in actual telecommunication sys- 
tems, particularly for mobile communication systems, has yet been suggested. The problems of applying conventional 
antenna array systems, which is based on the method of Eigen-Decomposition, is mainly due to its complexity and 
operating speed which is too large for real-time processing in telecommunication siystems. 

The conventional technique about the antenna array system was introduced in the following references: 

[1] M. Kaveh and A. J. Barabell, "The Statistical Performance of the MUSIC and Minimum-Norm Algorithms for 
Resolving Plane Waves in Noise," IEEE Trans., Acoust., speech and signal process., vol. ASSP-34, pp. 331-341 , 
April 1986. 

[2] T. Denidni and G. Y. Delisle, "A Nonlinear Algorithm for Output Power Maximization of an Indoor Adaptive 
Phased Array," IEEE Electromagnetic Compatibility, vol. 37, no. 2, pp. 201 -209, May, 1995. 

[3] B. G. Agee, S . V. Schell, and W. A. Gardner, "Spectral Self-Coherence Restoral: A New Approach to Blind Adap- 
tive Signal Extraction Using Antenna Arrays", Proc. of IEEE, Vol. 78, No. 4, pp. 753-767, Apr 1990. 

The problems in most conventional methods of designing antenna array system are, first, it, (except the method 
introduced in [3]), require some knowledge about the location of the wanted signal apriori, and second, it requires so 
many computations that the real-time processing cannot be performed. Especially, when the arrival angle of the wanted 
signal or the total number of signal sources is unknown, the required amount of computation becomes even larger, 
which makes it impossible to apply the conventional methdd of synthesizing the antenna array system to a practical sig- 
nal environment, such as mobile communications. Another undesirable feature of most conventional methods of 
designing antenna array systems is that the performance and/or the complexity of the system to be built is affected by 
the coherence and/or cross correlation of the wanted signal with respect to the interfering signals. This means that the 
antenna array system often requires additional complexities when the signals are fully or partially coherent so that the 
resultant system becomes too complicated for real-time processing, which is very important, especially in mobile com- 
munications. Another conventional methods such as the one shown in [3], does not work at all if the wanted signal can- 
not be extracted from the interfering signals with a particular frequency separation at a proper time lag. 

This invention introduces a new signal processing technology of designing an antenna array system that provides 
for a nice beam pattern having its maximum gain along the direction of the wanted signal maintaining the gain along the 
direction of interfering signals in a relatively much lower level: Under an assumption that the wanted signal is sufficiently 
larger in magnitude than each interfering signals, the proposed technique generates the desired beam pattern without 
requiring any knowledge about the wanted signal as well as the interfering signals. 

Another important and attractive aspect of the proposed technique is that the total required amount of computation 
is so small that the optimal parameters of the antenna array system are produced on a real-time basis. In fact, the signal 
processing apparatus, which forms the beamforming module of the antenna array system introduced in this invention, 
can easily be implemented with a normal, off the shelf digital signal processor. 

The primary objective of this invention is to introduce a new method of designing a signal processing apparatus, 
i.e., the beamforming module of an antenna array system, in order to apply it at the base station of a mobile communi- 
cation system.for receiving and transmitting the. signal of each subscriber in a cell with. a nice beam pattern which is 
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provided individually for each subscriber of the cell. The proposed technique can also be applied in other signal envi- 
ronments such as WLL(wireless local loop) and other fixed communications as well as mobile communications. 

SUMMARY OF THE INVENTION 

5 

It is an object of the present invention to provide a signal processing apparatus and method for enhancing the com- 
munication quality and increasing the communication capacity by reducing the effects of interference and noises with 
the nice beam pattern. 

And, the inventive signal processing apparatus and method introduce a simplified computational technique for gen- 
70 erating the nice beam pattern having its maximum gain along the direction of the wanted signal and maintaining the 
gain toward the direction of the interfering signals in as low a level as possible.. 

. To accomplish the object of the present invention, there is disclosed a signal processing apparatus for minimizing 
interference and for reducing effects of noise by controlling beam patterns of a telecommunication system having an 
antenna array, comprising: a means for computing. a parameter, gamma (y(/c)), by utilizing a predetermined adaptive 
15 gain (u), a signal vector (x(t)), each element of which is obtained from received signals at a corresponding antenna ele- 
ment and a final array output signal (y(t)) at the present snapshot; and a means for updating a gain vector (w) by utilizing 
said gamma {y(k)), the present value of said gain vector (w), said adaptive gain (n), said signal vector (x(t)) and said 
final array output (y(t)). 

In another aspect of the present invention, there is disclosed a- signal processing method for minimizing interfer- 
20 ence and reducing effects of noise by controlling beam patterns of a telecommunication system having an array 
antenna, comprising the steps of : (a) computing gamma (y) by utilizing said adaptive gain (^i), said signal vector (x) and 
said final array output signal (y) at the present snapshot; and (b) updating gain vector (w) by utilizing said gamma, the 
present value of said gain vector, said adaptive gain, said signal vector and said final array output. 

In another aspect of the present invention, there is disclosed a signal processing apparatus for minimizing interfer- 
25 ence and for reducing effects of noise by controlling beam patterns of a telecommunication system having an array, 
antenna, comprising: a means for generating an autocorrelation matrix (R) of received signals by utilizing said signal 
vector (x) at every snapshot; a means for computing said gamma (y) by utilizing said adaptive gain (u), the present value 
of said gain vector (w) and said autocorrelation matrix (R) at each snapshot; and a means for updating said gain vector 
(w) by utilizing said gamma (y), the present value of said gain vector (w), said adaptive gain (n) and said autocorrelation 
30 matrix (R). 

Also, in another aspect of the present invention, there is disclosed a signal processing method for minimizing inter- 
ference and for reducing effects of noise by controlling beam patterns of a telecommunication system having an array 
antenna, comprising the steps of: (a) generating an autocorrelation matrix (R) of received signals by utilizing said signal 
vector (x) at every snapshot; (b) computing gamma (y) by utilizing said adaptive gain fa), the present value of gain vec- 
35 tor (w) and autocorrelation matrix (R) at each snapshot; and (c) updating said gain vector (w) by utilizing said gamma 
(y), the present value of said gain vector (w), said adaptive gain fa) and said autocorrelation matrix (R). 

BRIEF DESCRIPTION OF THE DRAWINGS 

40 the novel features believed characteristic of the invention, as well as other features and advantages thereof, will 
best be understood by reference to the following detailed description of a particular embodiment, read in connection 
with the accompanying drawings, wherein: 

FIG. 1 is a block diagram of the signal processing apparatus according to the f irst embodiment of the present inven- 
ts tion. 

FIG. 2 is an example of the specified structure of the gamma-computing part shown in Fig. 1 ; 
FIG. 3A is an example of the specified structure of the gain vector updating part shown in Fig. 1 ; 
FIG. 3B is an another example of the specified structure of the gain vector updating part shown in Fig. 1 ; 
FIG. 4 is a block diagram of the signal processing apparatus according to the second embodiment of the present 
so invention. 

FIG. 5 is an example of the specified structure of the gamma-computing part shown in Fig. 4; 
FIG. 6A is a functional block diagram of the gain vector updating part shown in Fig. 4; 
FIG. 6B is another functional block diagram of the gain vector updating part shown in Fig. 4; and 
FIG. 17 shows a schematic block diagram of a telecommunication system that utilizes the signal processing appa- 
55 ratus according to the present invention shown in Fig. 1 or 4. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

A preferred embodiment of the present invention will be explained below with reference to the accompanying draw- 
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ings. 

The signal processing apparatus that is proposed in this invention generates a beam pattern having its maximum 
gain along the direction of the wanted signal maintaining the gain to the other directions in as low a level as possible. 
This can be accomplished by one of two approaches. 

5 The first approach is to optimize the value of the complex gain that is to be multiplied to each signal received at 

each antenna element, and the other approach is to optimize the value of the phase delay that is to be added to each 
signal received at each antenna element. Since each element of the gain vector in the first approach is to be weighted 
(multiplied) to each element of the signal vector, the gain vector is often referred to as the "weight vector" as well. 
Only the first approach is explained in detail in this document because of the following two reasons: first, these two 

10 approaches are mathematically equivalent, second, the second approach is more costy in hardware and the beamform- 
ing performance of the second approach is inferior to the first one in terms of angle accuracy. However, the second 
approach can easily be implemented by utilizing the procedure explained in this invention. Therefore a signal process- 
ing apparatus or signal processing method that is implemented in the second approach should be a part of this inven- 
tion as long as it utilizes the idea introduced in this document. 

75 In other words, this invention determines the complex gain vector "w" in such a way that the desired beam pattern 

be formed, and as a result the output of the array antenna system, i.e., the Euclidean inner product of the signals 
induced atthe antenna elements and the complex gain vector, should be as close to the wanted value as possible. 

If the magnitude of every element of the complex gain vector is normalized to 1 , multiplying the signal received at 
each antenna element by the corresponding element of the complex gain vector w is equivalent to adding the phase 

20 delay to the signal by the amount of the phase term of each corresponding element of the complex gain vector. There- 
fore, multiplying the signal vector by the gain vector is equivalent to adding the phase of the signal vector by the amount 
of the phase term of the gain vector. 

As mentioned above, the same effect can also be obtained by appending the time delay to the signal received at 
the i_th antenna element by the amount of <|>j divided by 2 n f c . where <h and f c denote the phase delay to be added 

25 to the signal received at the i_th antenna element and the carrier frequency, respectively. 
For a linear array system having a uniform spacing of 

2 

30 

between adjacent antenna elements, where X c denotes the wavelength at the carrier frequency, the signal induced at 
the m„th antenna element can be represented after the frequency down conversion as follows: 

35 M 

j(m-1)7t sin G k _ 1X /H v 

X m (t)=XS k (t)e K k + n m (t) (1) 

k=1 

40 where e k denotes the incident angle of the.kjh signal and S k (t) is the kjh transmitted signal observed at the 

receiving end. 

The subscript m in equation (1) represents the antenna element. The reference antenna element is assigned to be 

m=1 and the other antenna elements are assigned the next numbers, i.e., m=2, 3, in the order of the magnitude of 

the phase of the signal induced at each antenna element. 
45 In eq. (1), one of the M signals is the wanted signal. For example, when the S-,(t) is the wanted signal, the S-,(t) 

must be received at the antenna array system while all the other M-1 signals, i.e., S 2 (t), S 3 (t), S M (t), are interfering 

signals to be rejected together with the noise n m (t) for a good signal reception. 

Although eq. (1) is valid for the linear array with the uniform half -wavelength spacing, the technique provided in this 

invention can be generally applied to non-uniform spacing or non-linear array systems as well. 
so For non-uniform spacing arrays, if the distance of the m_th antenna element from the reference antenna element 

is d m , then there exists a phase difference in the signal induced at the m_th antenna element by 

2k ^ sine k 

55 

compared to the phase of the signal at the reference antenna element. Thus, the signal induced at the m_th antenna 
element for non-uniform and/or non-linear array systems can be written as follows: 
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M ]2nT JD sin 0 k 

X m (t) = £S k (t)e ° +n mW- 

k=1 



In this invention, in order to make the phase delay to be appended to each antenna element be a positive quantity, 
the reference antenna element is defined as the antenna element at which the induced signal has the latest phase in 
the receiving array. In the transmitting array system, therefore, the antenna element at which the induced signal has the 
earliest phase is the reference antenna element. 

w Defining the reference antenna element in the way explained above, the array antenna system can easily be 
designed by appending the zero phase delay to the signal at the reference antenna element and the proper positive 
amount of the phase delay to the signal at the other antenna elements. 

For an array antenna system consisting of N antenna elements, the array receives the N-by-1 signal vector at every 
snapshot. The autocorrelation matrix of the received signals can be written as shown in eq. (2). 

is The term "snapshot" in this document denotes the time period during which the new gain vector (or, phase delay 
vector) is computed upon receiving the new signal vector. In this invention, the array antenna system that adapts to the 
new signal vector can be designed at each snapshot by determining the proper gain vector (or, phase delay vector) for 
each new signal vector received at every snapshot. 



20 



40 



?x(J) = J Z?< t+1T s)? H < t+ 1T s ) - (2) 



25 where the underlined quantities denote the vector or matrix, T s is the snapshot period and superscript H is the 
Hermitian operator. The N-by-1 signal vector x(t), of which the number of elements is N, consists of the received signal 
x m (t)for m = 1, 2 N, which is explained in eq.(1) as follows: 

x(t)=[ Xl (t)x 2 (t). • -x N (t)] T . < 3 V 

30 

where superscript T denotes the transpose operator. 
However, eq. (2) is valid only when the arrival angles of all the signal components remain unchanged. In a time- 
varying environment where each signal source moves during the communication, as in the mobile communication envi- 
ronment, the autocorrelation matrix cannot be obtained by eq. (2) because the arrival angle of the signal source 
35 changes at every snapshot. 

Therefore, in time-varying signal environments, it is recommended that the autocorrelation matrix be computed in 
an iterative manner as follows: 



R x (J+1)=f-R x (J) + x((J+1)T s )x H ((J + 1)T s ) (4) 



where fix(J+1) and Bx(J) denote the autocorrelation matrix at the J+1st and J_th snapshot, respectively, and f 
denotes the forgetting factor in the range between 0 and 1 . 

Since communication environments, especially mobile communications, are generally time-varying environments, 
the autocorrelation matrix in this invention is computed by eq. (4) rather than eq. (2). 
45 From various computer simulations, it is recommended to set the value for the forgetting factor, 1, in the range 
between 0.8 and 0.99 for optimal performances in land mobile communications. 

Hereinafter, the design of the optimal array antenna system will be explained in more detail by taking the practical 
examples of the actual applications. 

The eigenvalues {A,/} of the autocorrelation matrix, determined by eq. (2) or (4), can be sorted by the magnitude as 
so X 1 > X 2 > * • * ^ 5in- Tne latest eigenvalue X.-, is determined by the signal components, not the noise components, 
regardless of the number of signal sources or antenna elements. 

Therefore, the eigenvector corresponding to the largest eigenvalue Xi exists in the signal subspace as follows: 

M 

55 e = X 4 i a(9 i) (5) 

where the complex quantity £ s is a constant determined by the magnitudes and distribution of the wanted and 
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interfering signals, and the vector a(9i) is the steering vector of the i_th signal component in the following form: 

m « jt sine , j(N-l)rc sin 9 = 

a(0i) = [1e • • - e ']. (6) 

Now, suppose the magnitude of the wanted signal is sufficiently larger than each of the interfering signals such that 
the condition shown in (7) is satisfied. 



I Si(t> I > > | Si(t) | for i*l 



(7 ) 



In a signal environment in which condition (7) is satisfied, the eigenvector ^ corresponding to the largest eigen- 
15 value can be approximated as: 

!i =A 1). (8) 

This means that the steering vector, a( 6 -,), of the wanted signal is almost the same as the eigenvector correspond- 
20 ing to the largest eigenvalue except that the complex-valued constant, % v is multiplied. 

Therefore, under the condition that the wanted signal is sufficiently larger than each of interfering signals, the max- 
imum gain of the array antenna system will approximately point to the direction of the source of the wanted signal if the 
gain vector to be appended to the antenna elements of the array system is determined by the eigenvector correspond- 
ing to the largest eigenvalue of the autocorrelation matrix of the signals impinging upon the array system. 
25 In conclusion of the above discussions, this invention suggests that the gain vector can be determined by the fol- 

lowing equation: 



Now, the practical way of computing the optimal weight vector will be presented. 

As mentioned previously, under a particular signal environment where the wanted signal is sufficiently larger than 
each of interfering signals, an. array antenna system having the desired beam pattern, which provides the maximum 
gain along the direction of the wanted signal source, can be obtained by taking the weight vector w with the normalized 
eigenvector e-, corresponding to the largest eigenvalue ^ of the autocorrelation matrix. 

However, to obtain the autocorrelation matrix itself requires a lot of computations, as shown in eqs. (2) and (4). 
Moreover, it is not a simple task to compute the eigenvector corresponding to the largest eigenvalue of the matrix What 
makes the problem even more complicated is that the arrival angle of each signal changes at every snapshot in mobile 
communications such that the eigenvector to be obtained varies at every snapshot. 

Considering the above-mentioned difficulties, this invention introduces a method of computing the weight vector w 
with" the approximated value for the eigenvector e 1 in an iterative way. This means that the weight vector w is computed 
by updating the solution of the previous snapshot through the iterative means as follows: 

w(k+1>w(k) +P (k)v(k) . (10) 

where the independent variable k is the time index representing the snapshot number, and p (k) and v(k) are the 
adaptive gain and search direction vector, respectively. Note that the gain vector w(k+1 ) shown in equation (1 0) should 
be normalized at each snapshot to make the magnitude of the gain vector.be 1 . 
so From equation (10). it is observable that the solution at the present snapshot can be obtained by updating the solu- 

tion of the previous snapshot in the direction indicated by v(k) by the amount indicated by p (k). 

In order to compute the solution for the gain vector in the iterative manner mentioned above, however, the answers 
for the following two questions must be satisfied: 

ss First, how do we set the initial value of the gain vector w(0) in the beginning? 

Second, how do we set the adaptive gain p (k) and the search direction vector y(k) at each snapshot? 

In this invention, the initial value of the gain vector w(0) is determined from the received signal vector x(0) as fol- 
lows: 1 
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5 where x-|{0), i.e., the first element of the signal vector x(0), is the signal induced at the reference antenna ele- 

ment at the very first snapshot. 

The reason why the vector j/v(0) can be determined by the equation (11) is that the received signal vector itself x(0) 
must be a good approximation for the searching eigenvector because the rank of the matrix at the initial snapshot is 1 , 
such that the number of the distinct nonzero eigenvalue is only 1, which must correspond to the signal received at the 
10 very first snapshot if the signal to noise ratio (SNR) is reasonably high. 

It is recommended that the magnitude of the adaptive gain does not exceed the reciprocal of N times of the aver- 
age power of the input signals, in order for the entire procedure of designing the antenna array system to converge, 
where N denotes the number of antenna elements. 

Starting from the initial value shown in (11), this invention introduces a new technique of designing the antenna 
is array system by updating the weight vector in the manner shown in equation (10). As mentioned previously, the key 
parts in updating the weight vector, as shown in equation (10), is to determine the search direction vector y(k) and the 
adaptive gain p (k). 

In order to find the search direction vector, let's consider a cost function defined as follows: 
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50 



f(w)= w(k) H R x w(k) + y(1-w(k) H w(k)) with a constraint | w(k)| 2 =1. (12) 



As can be verified mathematically, the maximum eigenvalue of the autocorrelation matrix R x can be.obtained by 
finding the maximum value of the functional (12), and the vector w corresponding to the maximum value of the func- 
25 tional (12) is an eigenvector corresponding to the maximum eigenvalue of the matrix R^ 

Since the gain vector w of the antenna array system should be determined by the eigenvector corresponding to the 
maximum eigenvalue in order to form a nice beam pattern having its maximum gain along the direction of the target sig- 
nal source, a search direction vector that maximizes the cost function (12) is to be found. 

The desired search direction vector described above can be obtained by setting the gradient of the function . (12) 
30 with respect to the weight vector to be zero as follows: 

V = 2 R x w - 2yw = 0. (13) 

In order to compute the gain vector satisfying equation (13) in an iterative manner as shown in equation (10), the 
35 search direction vector v{k) is set utilizing the result of equation 1(13) as follows: 

w(k+1) = w(k) + gp. v = w{k) + |i v{k) (14) 
= w(k) + »IR- X {k)wlk j - y{k)w{k)] 
= [(1-M*))/+^^W- 



It can be observed from equation (14) that the value for gamma y should be computed at each snapshot in order to 
45 obtain the gain vector w. The optimal value for gamma can be computed by substituting equation (14) into the constraint 
of equation (12) to result in the following expression for gamma: 



y 2 {k) , 2 [- + w H (k) R x {k)w[k)] y(k) + w H (k) R 2 x {k)w{k)+ ? w H (k) R x (k)w(k) = 0. (15) 



Putting two terms of equation (15) as follows 

a = [±+w H {k)R x {k)w(k)l 

55 fi ' 

b = w H {k) R z x (k)w(k) + ? w H (k) R x (k)w(k). 
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the value for gamma satisfying equation (15) can be obtained as 

y(/c)= a ± Ja^b. 

5 The smaller value for gamma is selected as follows: 

y(k) = a-Ja^b. (16) 

The entire procedure of computing the gain vector can be summarized as follows: 
w At the initial stage, the gain vector is set with the signal vector at the very first snapshot as w(0) = x(0)/x 1 (0) . 

The initial autocorrelation matrix is also set with the initial signal vector as R x (0) = x(0) x (0) . 

As the snapshot continues, the autocorrelation matrix is updated upon the reception of a new signal vector based 
on equation (4). The Value for gamma and gain vector are obtained according to equations (16) and (14), respectively. 
The update of the autocorrelation matrix, garnma, and gain vector is repreated with the new signal vector at every snap- 
15 shot. - 

According to the technique presented in this invention, the entire procedure of computing the gain vector and 
obtaining the final array output with the computed gain vector at each snapshot is tremendously simplified. The simpli- 
fication of the proposed method is mainly due to the fact that the technique disclosed in this invention does not require 
any apriori information regarding the location of the target signal source or interfering signal sources. Consequently, the 
20 proposed technique makes it possible to perform real-time processing of reception and transmission of the signals in 
most practical signal environments, such as mobile communications, by utilizing an ordinary, off the shelf digital signal 
processor (DSP). 

In fact, as shown in equations (4), (13) and (15), the total computational load for obtaining the gain vector is about 
0(2N 2 +6N) at each snapshot. It has been conf irmed in various computer simulations that the computational load of 
25 0(2N 2 +6N) is small enough to perform real-time processing of computing the gain vector and finally generating the 
array output utilizing a general -purpose DSP as long as the relative speed of the transmitter and receiver does not 
exceed 150km/h, as in land mobile communication. 

Although the amount of computation of the proposed technique has been reduced tremendously compared to the 
conventional methods, the complexity of this procedure increases in the second order of N, as the number N of antenna 
30 elements increases. The main reason for this is that the procedure includes the matrix operations. 

This invention discloses another technique of reducing the required amount of computation by setting the forgetting 
factor in computing the autocorrelation matrix with a particular value. 

Suppose the forgetting factor in equation (4) is set to be 0. This particularly means that the autocorrelation matrix 
is computed only with the instantaneous signal vector at each snapshot. Then, the entire procedure presented previ- 
35 ously is much more simplified. 

As a matter of fact, when the incident angle of the input signal varies significantly at each snapshot, it is not useful 
to consider the signal vector of the previous snapshots in computing the autocorrelation matrix. 

First, the computation of the autocorrelation matrix is simplified as follows: 

40 R x (J)«x(J)x H (J). (19) 

Substituting equation (19) into equations (1 3) through (1 5), the update of the search direction vector, gamma and 
gain vector is also simplified, respectively, as follows: 

45 v(k) = y*{k) x{k)-y(k) vv(/c), (20) 

y(7c) = a - Ja 2 -b , 

where, 

50 

a = 1 + |y(/c)| 2 , b = |y(/c)| 2 [||x(/c)|| 2 + ?]. 

and 

55 ' w</c+1) =. ^(/c) + v(/f) = + ^y*(/f) (21) 

Note that y(k) is the array output at the k_th snapshot, which is defined as y(k) = w H (k)x{k) . 
As shown in equation (21), if the forgetting factor is fixed at zero, then, since the matrix is determined by the signal 
vector of the present snapshot only, the procedure of computing the optimal weight vector is considerably simplified. 
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Moreover, the computation of the matrix at each snapshot is not needed at all, which means the calculation of equation 
(4) vanishes out of the entire procedure. 

From the numerical results obtained in the computer simulations, the proposed method, which accounts for the last 
previous signal vectors as well as for the present signal vector lor computing the autocorrelation matrix at each snap- 
shot, provides about a 12-t5dB improvement in SIR (signal-to-interference ratio), whereas the noise power is reduced 
by the number of antenna elements , i. e., the SNR (signal-to-noise ratio) is increased by the factor of N. 

On the other hand, the other method, which uses only the instantaneous signal vector at each snapshot, provides 
almost the same amount of improvement according to the noises, while about a 10-1 2dB improvement is obtained in 
terms of the SIR (signal-to-interference ratio). 

Consequently, the simplified version of the proposed method, which uses the signal vector at the present snapshot 
only, causes a degradation in SIR performance by about 2-3dB compared to the original version of the proposed 
method which uses the signal vectors of the previous snapshots as well as the current signal vector in computing the 
autocorrelation matrix. However, since the complexity of the entire procedure is tremendously reduced, a simplified ver- 
sion would cause a much easier implementation and cost reduction. 

Designing the array antenna system utilizing a simplified method, all the operations requiring the computational 
load of 0(N 2 ) disappear and the total computational load of the entire procedure becomes about 7N+16 additions and 
4N+7 multiplications. 

In order to implement the total system, which encounters both receiving and transmitting modes, the optimal weight 
vector computed during the receiving mode can be applied to obtain the optimal parameters for the transmitting mode. 

As mentioned previously, when the proposed signal processing apparatus, which provides the desired beam pat- 
tern, is adopted at the cell-site antenna system, we can achieve not only an increase of the channel capacity and an 
enhancement of the communication quality but also a considerable extension of the battery's life with each subscriber 
in the cell. 

An extension of the battery's life with each subscriber can be achieved because the cell-site antenna system adopt- 
ing the proposed beamforming technique provides much better communication efficiency by generating the main lobe 
along the direction of the wanted signal source compared to the conventional cell-site antenna system. 

Therefore, it is possible to perform an acceptable communication even with much less transmitting power at each 
subscriber's end. To reduce the transmitting power at each subscriber directly causes the life extension of the battery 
at each of the subscribers. 

Now, an explanation of the proposed apparatus and method in more detail will follow by taking practical examples. 
Embodied Example 1 

In this embodied example, a signal processing apparatus is introduced which computes the gain vector in real-time 
in order to generate the optimal beam pattern at the telecommunication system that employs the array antenna system. 
In this example, the autocorrelation matrix is updated with the instantaneous signal vector at each snapshot based on 
(19). Therefore, the autocorrelation matrix is actually not computed. As mentioned previously, the. gain vector is 
obtained from equation (21). The implementation of the signal processing apparatus is explained in detail as the first 
embodied example as follows. 

FIG. 1 is a block diagram of the signal processing apparatus according to an embodiment of the present invention. 

As illustrated in the figure 1, the signal processing apparatus according to the first embodiment of the present 
invention comprises a gamma computing part 1 1 for computing the gamma (y(/c)) and a gain vector updating part 12 for 
updating the gain vector (w). 

The gamma computing part 1 1 synthesizes the gamma (i(k)) by using a predetermined adaptive gain (u.), a signal 
vector (x(t)), each element of which is obtained from the received signal at the corresponding antenna element, and a 
final array output signal (y(t)) at the present snapshot. 

The gain vector updating part 12 updates the gain vector (w) by utilizing the gamma (y(/c)), the present value of the 
gain vector (w), the adaptive gain (n), the signal vector(x(t)) and the final array output (y(t)) at the present snapshot. 

Each part described above can be implemented in a form of software in a computing system, as well as in a form 
of hardware that is properly designed according to the procedures given in this invention. The ultimate goal of the signal 
processing apparatus is to generate the the gain vector (w) providing the optimal beam pattern for the telecommunica- 
tion system that employs the array antenna to produce the final array output sjgnal (y(t)) by computing the inner product 
between the signal vector received at the present snapshot and the gain vector (w). The details of computing the inner 
product is shown in Fig. 5. - . 

Fig. 2 illustrates an example of the specified structure of the gamma computing part 1 1 , which is a part of the signal 
processing apparatus shown in Fig. 1 . 

As shown in Fig. 2, the gamma computing part 1 1 comprises the following parts: a multiplying part G1 for cornput- 
ing the squared value of the magnitude of the final array output (y(t)); an adding part G2 for adding the result of said 
multiplying part G1 to the reciprocal ) of said adaptive gain; a multiplying part G3 for computing the squared value of 
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A, where A denotes the result of said adding part G2; a plurality of multiplying parts G4 for computing the squared value 
of the magnitude of each element of said signal vector (x); an adding part G5 for adding up alt the results of said multi- 
plying parts G4; an adding part G6 for adding the result of said adding part G5 to two-times ( j= ) of the reciprocal ( J ) of 
said adaptive gain (ji); a multiplying part G7 for multiplying the result of said adding part G6 by the result (|y| 2 ) of said 

5 multiplying part G1 ; an adding part G8 for subtracting B from the result A 2 of said multiplying part G3 where B denotes 
. the result of said multiplying part G7; a part G10 of computing a square root for generating the square root of the result 
of said adding part G8; and an adding part G9 for subtracting the result of said square-root-computing part G10 from 
the result of said adding part G2. 

The above-described parts can be implemented by means of software in a computing system or by means of hard- 

10 ware. 

The entire procedure can be summarized on a step-by-step basis as follows. 

The squared value of the magnitude of said final array output is obtained (first step). The reciprocal (J ) of said 
adaptive gain (u) is added to the result of G1 (second step). The squared value of A is computed, where A denotes the 
result of G2 (third step). The squared value of the magnitude of each element of said signal vector (x) is computed 

is (fourth step). All the results of G4 are added up (fifth step). The result of G5 is added to two-times the reciprocal of said 
adaptive gain (u). i.e. ( - ) (sixth step). The result of the procedure of G6 is multiplied by the result (|y| 2 ) of G1' (seventh 
step). B is subtracted from the result (A 2 ) of G3 where B denotes the result of G7 (eighth step). The square root of the 
result of G8 is computed (ninth step). And the result of G9 is subtracted from the result of G2, and thus, the resultant 
value of said gamma is produced as 

20 ■ - _ — : . 

y = A - JA Z -B 

(tenth step). 

Through the procedures presented above, the resultant gamma can also be viewed as 

25 

y = A - JA 2 -B 

where A and B denote 

30 - 1 H *'„ 

- + w R w 

M - - ~ 



and 

35 

!?" f? 2 - + \ - H 5 - ' 

40 respectively, with R being the autocorrelation matrix. 

Fig. 3a illustrates an example of the specified structure of the gain vector updating part 12, which is a part of said 
signal processing apparatus shown in Fig. 1 . 

As shown in Fig. 3a, the gain vector updating part comprises the following procedures: a multiplying part P1 for 
multiplying said gamma y by said adaptive gain \l ; an adding part P2 for subtracting the result of said multiplying pro- 
45 cedure P1 from 1 ; a plurality of multiplying parts P3 for multiplying the present value of each element of said gain vector 
by the result of said adding procedure P2; a multiplying part P4 for multiplying the complex conjugate of said final array 
output by said adaptive gain; a plurality of multiplying parts P5 for multiplying each element of said signal vector (x) by 
the result of said multiplying part P4; and a plurality of adding parts P6 for adding each output of said multiplying parts 
P3 to the corresponding output of said multiplying parts P5. 
so The outputs of the plural adding parts P6 form the final output of said gain vector updating part at each snapshot. 

The above-described procedures can also be implemented by means of software in a computing system as well as 
by means of hardware. 

The entire procedure can be summarized on a step-by-step basis as follows. 

The gamma (y) is multiplied by said adaptive gain (n) (first step). The result of said PI is subtracted from 1 (second 
55 step). The present value of each element of said gain vector is multiplied by the result of P2 (third step). The complex 
conjugate of the present value of said final array output (y) is multiplied by said adaptive gain (u) (fourth step). Each 
element of said signal vector at the present snapshot is multiplied by the result of P4 (fifth step). And each result of said 
P3 is added to each result of said P5, and thus, said gain vector is updated by w <~ (1 -^y) w + \iy*x (sixth step). 
Fig. 3b illustrates another example of the specified structure of the gain vector updating part 12, which is a part of 
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said signal processing apparatus shown in Fig. 1 . 

The gain vector updating part shown in Fig. 3b includes procedures for normalizing the resultant gain vector in 
addition to all the procedures that were in said gain vector updating part shown in Fig 3a. 

Since the procedures P1 to P6 are exactly the same as those in Fig. 3a, only the additional procedures P7, P8 ( P9 
and P10 for normalization of said gain vector are explained. 

As shown in Fig 3b, the gain vector updating part further comprises the following parts: a plurality of multiplying 
parts P7 for computing the squared value of the magnitude of each output of said adding parts P6; an adding part P8 
for adding up all the outputs of said multiplying parts P7; a part of computing square root P9 for obtaining the square 
root of the result of said P8; and a plurality of dividing parts P10 for dividing each output of said P6 by the result of said 
P9. 

- The above-described parts can also be implemented by means of software in a computing system or by means of 
hardware. 

The entire procedure of Fig. 3b is for updating the gain vector by w <- (1-ny) w+ uy*x and then normalizing the 
resultant vector by 

w 
||w[| 

Therefore, the magnitude of the gain vector produced in accordance with Fig. 3b is always normalized to be 1 . 
Embodied Example 2 

Here, another example of implementing a signal processing apparatus is introduced which computes the gain vec- 
tor in real-time in order to generate the optimal beam pattern in the telecommunication system that employs the array 
antenna system. In this example, the autocorrelation matrix is computed with the signal vectors based on equation (4). 
By adopting a proper value for the forgetting factor, the autocorrelation matrix is updated at each snapshot. As men- 
tioned previously, the gain vector is obtained from equation (14). 

Although the complexity of the signal processing technique disclosed in this example is a little heavier than that of 
the first example, the performance is a little better in terms of a SIR or BER improvement. Therefore, the first example 
can be applied to a relatively simple system, whereas the second example can be applied to a larger system that 
requires more accuracy. The implementation of the signal processing apparatus is explained in detail as the second 
embodied example as follows. 

FIG. 4 is a block diagram of the signal processing apparatus according to an embodiment of the present invention. 

The signal processing apparatus according to the second embodiment of the present invention comprises an auto- 
correlation matrix updating part 20, a gamma computing part 21 and a gain vector updating part 22. 

The autocorrelation matrix updating part 20 computes a new value for the autocorrelation matrix at each snapshot 
upon the reception of a new signal vector (x), each element of which is obtained from the signal induced at the corre- 
sponding antenna element of the array system at every snapshot. 

The gamma computing part 21 computes the value for gamma(y) by utilizing said autocorrelation matrix(£7), said 
adaptive gain(u) and the present value of said gain vector(w). 

The gain vector updating part 22 updates the gain vector by utilizing said gamma(y), the autocorrelation matrixO, 
the adaptive gain(^) and the present value of said gain vector(w). 

As mentioned previously, every part included in the signal processing apparatus shown in Fig. 4 can be imple- 
mented by means of software in a computing system as well as hardware. The details in each part of the signal 
processing apparatus are shown in Figs. 5, 6a and 6b. 

FIG. 5 is a functional block diagram of the gamma computing part, illustrated in Fig. 4, according to an embodiment 
of the present invention. 

The gamma computing part according to the second embodiment of the present invention computes the value for 
said gamma by utilizing the autocorrelation matrix, the present value of the gain vector and the adaptive gain through 
a computational procedure described as follows. 

Each row of the autocorrelation matrix is multiplied by the present value of said signal vector to produce a vector 
E, i.e., E = R w (first step, 51) . The complex conjugate of each element of the gain vector is multiplied by the corre- 
sponding element of E, i.e., X = w H E (second step. 52) . The L 2 -norm of the vector E is computed, i.e., 
F = E H E (third step, 53) . Two-times of the reciprocal of the adaptive gain is multiplied by X. i.e., 



G = - k (fourth step, 54). 
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The reciprocal of the adaptive gain is added to X, i.e., 

A = ^ + X (fifth step, 55). 

5 

The squared value of A is computed, i.e., C = A 2 (sixth step, 56) . F is added to G, i.e., 
B = F + G (seventh step, 57) . And, the gamma is computed by y = A - VC^B (eighth step, 58) . 

FIG. 6a is a functional block diagram of the gain vector updating part, illustrated in Fig. 4, according to an embod- 
w iment of the present invention. 

The gain vector updating part according to the second embodiment of the present invention updates the value for 
the gain vector by utilizing the autocorrelation matrix, the present value of the gain vector, and the adaptive gain through 
a computational procedure described as follows. 

The adaptive gain is multiplied by the gamma, i.e., p 1 = u,y (first step, 61). The result of the first step 61 is sub- 
15 traded from 1 , i.e., p 2 = 1 - p ., (second step, 62) . The result of the second step 62 is added to the main diagonal ele- 
ments of matrix n R to produce resultant matrix. Q, i.e., Q = p 2 /+ u R (third step, 63). And, the gain vector is 
updated by w <r~ Q w (fourth step, 64) '-. 

FIG. 6b is a functional block diagram of the gain vector updating part, illustrated in Fig. 4, according to an embod- 
iment of the present invention. 

20 The gain vector updating part shown in Fig. 6b includes normalization procedures in addition to what is included in 
the gain vector updating part shown in Fig. 6a. 

The gain vector updating part according to the second embodiment of the present invention updates the value for 
the gain vector by utilizing the autocorrelation matrix, the present value of the gain vector, and the adaptive gain through 
a computational procedure described as follows. 

25 The adaptive gain is multiplied by the gamma, i.e., p 1 = ny (first step, 61) . The result of the first step 61 is sub- 
tracted from 1 , i.e.. p 2 = 1 - p n (second step, 62) . The result of the second step 62 is added to the main diagonal ele- 
ments of matrix ^ f?to produce resultant matrix Q, i.e., Q = p 2 / + n R (third step, 63) . The Q matrix is multiplied by 
the present value of the gain vector, i.e., D~Qw (fourth step, 64) . The magnitude of each element of the vector D is 
squared and the results of the squaring are added up together, i.e., 

N 

P3= D D /| 2 ( ffln STEP, 65). 
/=1 

35 

And, the vector D is divided by p 3 to produce the updated gain vector, i.e., 

D 

w< == (sixth step, 66). 

40 V p 3 ' 

The gain vector produced in Fig 6b is normalized for the magnitude of the resultant gain vector to be 1 at each 
snapshot. 

FIG. 7 shows a schematic block diagram of a telecommunication system that utilizes the signal processing appa- 

45 ratus according to the present invention shown in Fig. 1 or 4. 

In Fig. 7, the reference number 1 denotes an array antenna, 7 a receiving apparatus, 8 an inner product computing 
apparatus (which is sometimes denoted as the part of generating the final array output), and 9 the signal processing 
apparatus according to the present invention, respectively. 

As illustrated in the figure,, the telecommunication system comprises the following parts: the array antenna 1 (or, 

so called simply, "array", "antenna array", or, "array of antenna elements"), composed of the plural antenna elements 1 1 [ 
each of which is arranged by a predetermined geometry, that supplies the signal induced at each antenna element to 
the corresponding port of the receiving apparatus 7; the signal receiving apparatus 7 that generates the signal vector 
(x(t)) from the signals induced at each antenna element of the antenna array 1 by utilizing the proper signal-receiving 
parts, such as filtering, frequency-down-conversion and demodulation (or, quasi-quadrature detection); the inner prod- 

55 uct computing apparatus 8 for generating the final array output ( y(t)) by computing the Euclidean inner product between 
the two complex-valued vectors, ( y(t) = w x(t) ), i.e., the signal veictor (x(t)) produced from the receiving part 7 and 
the gain vector (w) provided from the signal processing apparatus 9; and the signal processing apparatus 9 that com- 
putes the gain vector (w) by processing the signal vector (x(t)) together with the final array output (y(t)) obtained at the 
last previous snapshot for the inner product computing apparatus 8 to generate the final array output (y(t)) at the present 
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snapshot 

The telecommunication system consists of the receiving apparatus 7, the signal processing apparatus 9 and the 
inner product computing apparatus 8 for generating the final array output. The receiving apparatus generates the signal 
vector (x(t)) from the signals induced at the antenna elements 1 1 through a conventional signal reception part,, such as 
a frequency-down-conversion and demodulation (or quasi-quadrature detection). 

When the technique provided in this invention is applied in the CDMA (Code Division Multiple Access) system, the 
receiving apparatus 7 includes a cross-correlation part for cross^correlating the received signal with the code sequence 
assigned to the wanted signal source. The signal vector (x(t)) obtained from the receiving apparatus 7 is sent to the sig- 
nal processing apparatus 9 and .the inner product computing apparatus 8. 

The signal processing apparatus 9 produces the optimal gain vector (w), which is sometimes referred to as the 
"weight vector", from the signal vector (x(t)) at the present snapshot and the finai array output (y(t)) computed at the last 
previous snapshot. The optimal weight vector (w) is sent to the inner product computing apparatus for the final array 
output (y(t)) of the next snapshot to be computed as a result of the inner product of the signal vector (x(t)) and weight 
vector (w), i.e., y(t) = w H x(t) . 

The key part of the terecommunication system shown in Pig. 7 is the signal processing apparatus 9 producing the 
optimal weight vector (x(t)), which gives the array antenna system the optimal beam pattern having its maximum gain 
along the direction of the wanted signal source and small gain to the direction of the interfering signal sources. 

As mentioned earlier, the signal processing apparatus or signal processing technique provided in this invention 
gives the following advantages: first, the communication capacity is increased as much as the signal-to-interference 
ratio is increased, and second, the communication quality is enhanced as much as the signal-to-noise ratio and the sig- 
nal-to-interference ratio is increased. The best feature of the proposed technique in this invention is that the required 
amount of computation to achieve all the merits is extremely small so that the proposed technique can be easily imple- 
mented with the normal digital signal processor in real-time processing. 

Although the specific embodiments of the present invention have been disclosed and described, it is apparent that 
those who skilled in the art will appreciate that various modifications, additions and substitutions are possible, without 
departing from the scope and the spirit of the present invention as disclosed in the accompanying claims. Therefore, it 
should be understood that the present invention is not limited to the particular embodiment disclosed herein as the best 
mode contemplated for carrying out the present invention. 

Claims 

1 . A signal processing apparatus for minimizing interference and for reducing effects of noise by controlling beam pat- 
terns of a telecommunication system having an array antenna, comprising: 

a means for computing a parameter, gamma (y(/c)), by utilizing a predetermined adaptive gain (u.), a signal vec- 
tor (x(t)), each element of which is obtained from received signals at a corresponding antenna element and a 
final array output signal (y(t)) at the present snapshot; and 

a means for updating a gain vector (w) by utilizing said gamma (y(/c)), the present value of said gain vector (w), 
said adaptive gain (u), said signal vector (x(t)) and said final array output (y(t)). 

2. The signal processing apparatus according to claim 1. wherein said gain vector (w) is determined by a value of an 
eigenvector corresponding to the maximum eigenvalue of an autocorrelation matrix of the signals received at each 
antenna element of said array antenna. 

3. The signal processing apparatus according to claim 2, wherein said gain vector (w) is determined by multiplying a 
predetermined constant on each element of said eigenvector corresponding to said maximum eigenvalue of said 
autocorrelation matrix, in order to modify said gain vector without changing beam-pattern characteristics of said 
eigenvector of said maximum eigenvalue. 

4. The signal processing apparatus, according to claim 2, wherein said gain vector (w) is determined by normalizing 
said eigenvector corresponding to said maximum eigenvalue of said autocorrelation matrix, such that a magnitude 
of the normalized eigenvector becomes 1 and a beam-papern characteristics of said eigenvector of said maximum 
eigenvalue remains unchanged. 

5. The signal processing apparatus according to claim 2, wherein said autocorrelation matrix is computed by adding 
a first term and a second term, as shown in the equation given below: (in the equation, said first term is the auto- 
correlation matrix at the last previous snapshot multiplied by a forgetting factor, of which the magnitude is between 
0 and 1 , and said second term is a signal matrix computed with said signal vector (x(t)) obtained from each antenna 
element of said array antenna at the present snapshot.) 
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R x (J+1)rf-B x (J) + x((J+1)T s )x H ((J+^ 

where R X (J+1) and R X (J) denote said autocorrelation matrix at the J+1_st and J_th snapshots, respectively, 
f is said forgetting factor of which the magnitude lies between 0 and 1 , T s is a snapshot period and superscript H 
denotes a Hermitian operator. 

The signal processing apparatus according to claim 5, wherein said forgetting factor in forming said autocorrelation 
matrix is set to be 0 for simplifying the entire procedure. 

The signal processing apparatus according to claim 2, wherein said eigenvector corresponding to said maximum 
eigenvalue is computed by the procedures of: 

(a) determining said gain vector with a normalized value of the received signal vector (2<(0)), wherein the phase 
of each element of said signal vector is modified in such a way that the resultant array output is synchronized 
to the phase of the signal received at said reference antenna element during the first snapshot; and 

(b) updating said gain vector of the last previous snapshot in such a way that a predetermined cost function, 
f(w) = w H ft x vv+ y(1 - w H w) , is maximized satisfying a constraint | w(/c)| 2 =1 at each snapshot, and a gain 
value to be multiplied to said signal received at said reference antenna element at each snapshot is maintained 
to be a real quantity during the second snapshot and on. 

The signal processing apparatus according to claim 7, wherein said reference antenna element is determined by 
an antenna element of which the phase of said signal is the latest of all said antenna elements in said array antenna 
at the present snapshot. 

The signal processing apparatus according to claim 7, wherein said reference antenna element is determined by 
said antenna element of which the physical distance from a signal source to be communicated with at the present 
snapshot is farthest compared to the other antenna elements in said array antenna. 

30 10. The signal processing apparatus according to claim 1, wherein said means for computing said gamma (y) com- 
prises: 

a first multiplying means for computing the squared value of the magnitude of the final array output (y(t)); 
a first adding means for adding the result of the first multiplying means to the reciprocal ) of said adaptive 
35 gain; 

a second multiplying means for computing the squared value of A, where A denotes the result of the first add- 
ing means; 

a plurality of multiplying means G4 for computing the squared value of the magnitude of each element of said 
signal vector (x); 

40 a second adding means for adding up all the results of said multiplying means; 

a third adding means for adding the result of the second adding means to two-times ) the reciprocal ) of 
said adaptive gain (|i); 

a third multiplying means for multiplying the result of the third adding means by the result (|y| 2 ) of the first mul- 
tiplying means; 

45 a fourth adding means for subtracting B from the result A 2 of the second multiplying means, where B denotes 

the result of the third multiplying means; 

a means for computing a square root for generating the square root of the result of the fourth adding means; 
and 

a fifth adding means for subtracting the result of said square-root-computing means from the result of the first 
so adding means to generate the valuie of the gamma by 

y(k) = a - Ja 2 -b 

55 where 

a = l + |y(/c)| 2 , 
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b = \y(k)\ 2 [\\x(k)\\ 2 + ^). 



11. The signal processing apparatus according to claim 1, wherein said means for updating said gain vector (w) com- 
prises: 

a first multiplying means for multiplying said gamma y by said adaptive gain (ji); " 
a first adding means for subtracting the result of the first multiplying means from 1 ; 

a" plurality of multiplying means for multiplying the present value of each element of said gain vector by the 
result of the first adding means; 

a second multiplying means for multiplying the complex conjugate (y*) of said final array output by said adap- 
tive gain; 

a plurality of multiplying means for multiplying each element of said signal vector by the result of the second 
multiplying means; and 

a plurality of adding means for adding each output of said multiplying means to the corresponding output of 
said multiplying means to provide the final output of said gain vector updating means at each snapshot as 
w <r- (1-uy) w+ \iy*x- 

12. The signal processing apparatus according to claim 1 1 , wherein the means of updating said gain vector further 
comprises: 

a plurality of multiplying means for computing the squared value of the magnitude of each output of said adding 
means; 

an adding means for adding up all the outputs of said multiplying means; 
a means for computing square root of the result of said adding means; and 

a plurality of dividing means for dividing each output of said adding means by the result of said square-root 
computing means to generate a normalized value for the gain vector as 

w 

- IMI 

13. A signal processing method for minimizing interference and reducing effects of noise by controlling beam patterns 
of a telecommunication system having an array antenna, comprising the steps of : 

. (a) computing gamma (y) by utilizing said adaptive gain'(jj), said signal vector (x) and said final array output 
signal (y) at the present snapshot; and 

(b) updating gain vector (w) by utilizing said gamma, the present value of said gain vector/said adaptive gain, 
said signal vector and said final array output. 

14. The signal processing method according to claim 13, wherein the step of (a) computing said gamma (y) comprises 
the following substeps of: 

(a-1 ) computing the squared value of the magnitude of said final array output; 

(a-2) adding the reciprocal ( 2 ) of said adaptive gain (p) to the result of said substep (a-1); 

(a-3) computing the squaredvalue of the result (A) of said substep (a-2); 

(a-4) computing the squared value of the magnitude of each element of said signal vector (x); 

(a-5) adding up all the results of said substep (a-4); 

(a-6) adding the result of said substep (a-5) to the multiplication of the reciprocal of said adaptive gain (ji) by 
2. i-e.,(i); 

(a-7) multiplying the result of said substep (a-6) by the result (|y| 2 ) of said substep (a-1). 
(a-8) subtracting the result (B) of said substep (a-7) from the result (A 2 ) of said substep (a-3); 
(a-9) computing the square root of the result of said substep (a-8); and 

(a-10) subtracting the result of said substep (a-9) from the result of said substep (a-2), and thus, the resultant 
value of said gamma is produced as 



y = A - Ja 2 -B. 
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15. The signal processing method according to claim 13, wherein the step of (b) updating said gain vector comprises 
the following substeps of: 

(b-1) multiplying said gamma (y) by said adaptive gain (u); 
(b-2) subtracting the result of said substep (b-1) from 1 ; 

(b-3) multiplying the present value of each element of said gain vector by the result of said substep (b-2); 
(b-4) multiplying the complex conjugate of the present value of said final array output (y) by said adaptive gain 

(n); 

(b-5) multiplying each element of said signal vector at the present snapshot by the result of said substep (b-4); 
and 

(b-6) adding each result of said substep (b-3) to each result of said substep (b-5), and thus, said gain vector is 
updated by w <r~ (I 1 uy) w + \iy*x. 

16. The signal processing method according to claim 13, wherein the step of (b) updating said gain vector comprises 
15 the following substeps of: 

(b-1) multiplying said gamma (y) by said adaptive gain (u); 
(b-2) subtracting the result of said substep (b-1) from 1 ; 

(b-3) multiplying the present value of each element of said gain vector by the result of said substep (b-2); 
(b-4) multiplying the complex conjugate of the present value of said final array output (y) by said adaptive gain 

(n); 

(b-5) multiplying each element of said signal vector at the present snapshot by the result of said substep (b-4); 
(b-6) adding each result of said substep (b-3) to each result of said substep (b-5); 
(b-7) computing the squared value of the magnitude of each output of said substep (b-6); 
(b-8) adding up all the outputs of said squared values; 
(b-9) computing the square root of the result of substep (b-2); and 
(b-10) dividing each output of said substep (b-6) by the result of substep (b-9) to generate 

w 

w 4— " 

- IMI" 

17. A signal processing apparatus for minimizing interference and for reducing effects of noise by controlling beam pat- 
terns of a telecommunication system having an array antenna, comprising: 

35 - . 

a means for generating an autocorrelation matrix (JR) of received signals by utilizing said signal vector (x) at 
every snapshot; 

a means for computing said gamma (y) by utilizing said adaptive gain (u), the present value of said gain vector 
(w) and said autocorrelation matrix (B) at each snapshot; and 
46 a means for updating said gain vector (w) by utilizing said gamma (y), the present value of said gain vector (w), 

said adaptive gain (\x) and said autocorrelation matrix (R). 

18. The signal processing apparatus according to claim 17, wherein said gain vector (w) is determined by a value of an 
eigenvector corresponding to the maximum eigenvalue of an autocorrelation matrix of the signals received at each 

45 antenna element of said array antenna. 

19. The signal processing apparatus according to claim 18, wherein said gain vector (w) is determined by multiplying 
a predetermined constant to each element of said eigenvector corresponding to said maximum eigenvalue of said 
autocorrelation matrix, in order to modify said gain vector without changing beam-pattern characteristics of said 

so eigenvector of said maximum eigenvalue. 

The signal processing apparatus according to claim 18, wherein said gain vector (w) is determined by normalizing 
said eigenvector corresponding to said maximum eigenvalue of said autocorrelation matrix, such that the magni- 
tude of the normalized eigenvector becomes 1 and a beam-pattern characteristics of said eigenvector of said max- 
imum eigenvalue remains unchanged. 

21. The signal processing apparatus according to claim 18, wherein said autocorrelation matrix computing part (20) 
computes the autocorrelation matrix by adding a first term and a second term, as shown in the equation given 
below: (in the equation, said first term is the autocorrelation matrix at the last previous snapshot multiplied by a for- 
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getting factor of which the magnitude is between 0 and 1,. and said second term is a signal matrix computed with 
said signal vector (x(t)) obtained from each antenna element of said array antenna at the present snapshot.) 

R X (J+1) = f • R X (J) + x((J+1)T s ) x H ((J+1)T s ) 

where Fyj+I) and R^J) denote said autocorrelation matrix at the J+1_st and J_th snapshots, respectively, 
f is said forgetting factor of which the magnitude lies between 0 and 1 , T s is a snapshot period and superscript H 
denotes a Hermitian operator. 

22. The signal processing apparatus according to claim 18, wherein said eigenvector corresponding to said maximum 
eigenvalue is computed by the procedures of: 

(a) determining said gain vector with a normalized value of the received signal vector {£€))), wherein the phase 
of each element of said signal vector is modified in such a way that the resultant array output is synchronized 
to the phase of the signal received at said reference antenna element during the first snapshot; and 

(b) updating said gain vector of the last previous snapshot in such a way that a predetermined cost function, 
i(w) = w H R x w+ y(1 ■ w H w) , is maximized satisfying a constraint \ w(k)\ =1 at each snapshot, and a gain 
value tobe multiplied to said signal received at said reference antenna element at each snapshot is maintained 
to be a real quantity during the second snapshot and on. 

23. The signal processing apparatus according to claim 22, wherein said reference antenna element is determined by 
an antenna element of which the phase of said signal is the latest of all said antenna elements in said array antenna 
at the present snapshot. 

24. The signal processing apparatus according to claim 22, wherein said reference antenna element is determined by 
said antenna element of which the physical distance from a signal source to be communicated with at the present 
snapshot is farthest compared to the other antenna elements in said array antenna. 

25. The signal processing apparatus according to claim 17, wherein said gamma computing part comprises: , 

a first multiplying means for multiplying each row of the autocorrelation matrix by the present value of said sig- 
nal vector to produce a vector E, i.e., E = R w ; 

a second multiplying means for multiplying the complex conjugate of each element of the gain vector by the 

corresponding element of E, i.e., X = w E ; h 

a first computing means for computing the L 2 -norm of the vector E, i.e., F= E .E ; - 

a third multiplying means for multiplying two-times of the reciprocal of the adaptive gain by X, i.e., 

G= - x: 

a first adding means for adding the reciprocal of the adaptive gain to X, i.e., 

A = - + X; 

a second computing means for computing the squared value of A, i.e., C = A ; 
a second adding means for adding the F to the G, i.e., B = F + G ; a nd 
a third computing means for computing the gamma by y = A - VC-B . 

26. The signal processing apparatus according to claim 17, wherein said gain vector updating part comprises: 

a first multiplying means for multiplying the adaptive gain by said gamma, i.e., p A = u.y : 

a subtracting means for subtracting the result of said first multiplying means from 1 , i.e., p 2 = 1 - p 1 ; 

a adding means for adding the result of said subtracting means to the main diagonal elements of a matrix \iR 

to produce a matrix (2 i.e., 0=p 2 /+uf?, where /denotes the identity matrix; and 

a second multiplying means for multiplying the matrix Q by the present value of the gain vector w to produce 
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the updated gain vector by w<r- Qw. 



27. The signal processing apparatus according to claim 17, wherein said gain vector updating part comprises: 

a first multiplying means for multiplying the adaptive gain by the gamma, i.e., p 1 = p.y ; 
a subtracting means for subtracting the result of said first multiplying means from 1, i.e., p 2 = 1 - p 1 ; 
a means for adding the result of said subtracting means to the main diagonal elements of a matrix R to pro- 
duce resultant matrix G, i.e., Q = p 2 / + u. R , where /denotes the identity matrix; 

a second multiplying means for multiplying the Q matrix by the present value of the gain vector, i.e., D = Q w ; 
a means for computing the L 2 -norm of the vector D, i.e., 



28. A signal processing method for minimizing interference and for reducing effects of noise by controlling beam pat- 
terns of a telecommunication system having an array antenna, comprising the steps of: 

(a) generating an autocorrelation matrix (R) of received signals by utilizing said signal vector (x) at every snap- 
shot; - 

(b) computing gamma (y) by utilizing said adaptive gain (ji), the present value of gain vector (w) and autocor- 
relation matrix (R) at each snapshot; and 

(c) updating said gain vector (w) by utilizing said gamma (y), the present value of said gain vector (w), said 
adaptive gain (^i) and said autocorrelation matrix (R). 

29. The signal processing method according to claim 28, wherein the step of (b) computing said gamma comprises the 
substeps of: 

(b-1) multiplying each row of the autocorrelation matrix by the present value of said signal vector to produce a 
vector E, i.e., E = R w ; 

(b-2) multiplying the complex conjugate of each element of the gain vector by the corresponding element of £, 
i.e., X = w E ; - 

(b-3) computing the L 2 -norm of the vector E, i.e., F = E H E ; 

(b-4) multiplying two-times of the reciprocal of the adaptive gain by X, i.e., 




and 

a means for dividing the vector D by p 3 to produce the updated gain vector, i.e., 



D 





(b-5) 



adding the 



reciprocal of the adaptive gain to X, i.e., 




(b-6) computing the squared value of A, i.e., C = A 2 ; 
(b-7) adding the F to the G, i.e., B = F + G ; and 
(br8) computing the gamma by y = A - VC-B . 



30.. The signal processing apparatus according to claim 28, wherein the step of (c) updating said gain vector comprises 



the substeps of: 
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(c-1) multiplying the adaptive gain by said gamma, i.e., p 1 =>y ; 

(c-2) subtracting the result of said substep (c-1) from 1, i.e., p 2 = 1 " P i ■ 

(c-3) adding the result of said substep (c-2) to the main diagonal elements of a matrix \xR\o produce a matrix 
Q, i.e., Q = p 2 /+}ift where /denotes the identity matrix; and 
5 (c-4) multiplying "the matrix Q by the present value of the gain vector w to produce the updated gain vector by 

w <r- Q w. 

31 . The signal processing apparatus according to claim 28, wherein the step of (c) updating said gain vector comprises 
the substeps of: 

10 

— (c-1) multiplying the adaptive gain by the gamma, i.e., p 1 = jxy , . 

(c-2) subtracting the result of said substep (c-1) from 1, i.e., p 2 = 1 - p n ; 

(c-3) adding the result of said substep (c-2) to the main diagonal elements of a matrix p. f?to produce resultant 
matrix Q, i.e., Q = p 2 / + H R where /denotes the identity matrix; 
is (c-4) multiplying the Q matrix by the present value of the gain vector, i.e., D=Qw\ 

(c-5) computing the L 2 -norm of the vector D, i.e., 

p 3 = fjlA! 2 ; 

20 

and 

(c-6) dividing the vector £> by p 3 to produce the updated gain vector, i.e., 
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FIG. 2 
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FIG. 3B 
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FIG. 4 
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FIG. 5 
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FIG. 6B 
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FIG. 7 
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